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Abstract 

The pseudoscalar quarkonia exclusive decays to light mesons still poses a challenge to the theo- 
retical understanding of quarkonium properties in decay. In this work, we evaluate the processes of 

pseudoscalar heavy quarkonium decays into vector meson pairs, especially the helicity suppressed 
(N 
^ ' processes of r]b — )■ J/tpJ/ip and rjc — >• VV. In the frame of NRQCD, the branching fraction of 
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Br[r]h — >• J/tpJ/^] are evaluated at the next-to-leading order of perturbative QCD; and within 
the light-cone distribution formalism, we calculate also the higher twist effects in these processes. 

in 

^D I Numerical results show that the higher twist terms contribute more than what from the NLO QCD 
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corrections in the process of r/f, — )• J /ipJ/il). It is found that the experimental results on rjc — )• VV 
are hard to be understood by merely the quark model and perturbative QCD calculation. 
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I. INTRODUCTION 

In high energy physics, heavy quarkonium study is one of the most interesting fields 
and it plays an important role in the understanding of the configurations of hadrons and 
the nonperturbative behavior of strong interaction. On one hand, the heavy quark masses 
enable the perturbative QCD(pQCD) calculation for quarkonium production and decay 
possible. On the other hand, due to the non-relativistic nature of heavy quarkonium, one 
may investigate their properties through a more transparent way, i.e. the effective theory of 
non-relativistic QCD(NRQCD) ll|. 

It is well known that the S-wave spin-triplet heavy quarkonium states, the J /ip and 
T, can be produced directly in e^e~ annihilation, and be measured via lepton pair decay 
mode distinctively. These characters lead to rich experimental data and deep investigations 
on them. While for their spin-singlet partners, the rji, and rjc-, things are not that easy. 
At present, people know relatively much less about their properties, especially for rji,. For 
rjci though there have been many measurements in experiment on its various decay modes, 
puzzles remain in confronting theoretical explanations to the experimental data, such as in 

nn 

r]c decay to vector meson pair [2|, |3[. For r/^ study, there have many theoretical scenarios 



been put forward |4j-|l0|. and several experiments been conducted |llHl5|. However, it was 
fixed only in very recently by BaBar collaboration through T(35') — )■ T^t + 7 process 16| and 
later on confirmed by CLEO-c experiment jlTJ]. About 77^ so far we merely know the mass, 
its other properties are remaining unclear and waiting for further investigations. It is worth 
noting that both Babar and CLEO-c measurements are indirect ones. For further study on 
?7fe physics, direct measurements on its decay products are necessary. 

For the direct rjh detection, Braaten et al. suggested to measure its exclusive decay 
products, the J /if) pair in 77^ decays [6|]. In comparison with the experimental result for rj^. — >■ 
(f)(j) and by some scaling assumptions, they estimated the branching ratio of r/^ — )■ J /ipj /ip 
mode to be 7 X 10^^^^, which hence is promising to be observed in the Fermilab Tevatron 
Run II experiment. So far there has been no conclusive report from the experiment yet, and 
the theoretical estimation was questioned by Maltoni and Polosa [[18J . In the expectation of 



helicity conservation rules 19|], the leading order calculation in the nonrelativistic limit gives 



null result. The calculations on next-to-leading order QCD \J\ and relativistic corrections 



8| both yield the branching ratios to be about 10~*. Recently, Braguta et al. reevaluate 

the r^fc — 7- J /ipj /ip process in the light cone formalism and find that after including the next- 

;o-leading twist contribution the branching fraction can be as large as (6.2 ± 3.5) x 10^^ 
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9j. The authors of Ref. |9| claim that the result in [8| dose not agree with theirs. The 



form factor obtained in Ref. 



Is 



contains double logarithms, whereas they are absent in [9|. 
To carry on an independent calculation of the ?7b — > J /ipJ /ip process is therefore one of the 
aims of this work. 

Similarly, the processes rjc — )■ VV are also governed by the helicity selection rules, but 
experiment gives extremely large results [^ , which stands as a long term puzzle existing in 
the charmonium physics. The higher order radiative corrections give this issue no help, since 
they are all suppressed by the light quark masses. Although beyond the scope of pQCD, 
some nonperturbative models have been put forward and considered to be solutions to the 



problem, such as the intermediate meson exchange model 20|] and the charmonium light 
Fock component admixture model [21|], to further investigate it in pQCD is still necessary. 

The rest of this paper is organized as follows: in section II, we calculate the branching 
ratio of the process rj^ — )■ J /ipj /tp at one-loop level; in section III, we evaluate the higher twist 
effects in processes r^^ — )■ J jipj /ip and r/c — > VV] in section IV, summary and conclusions 
are presented. 

II. NLO QCD RESULT FOR 7?^ -^ J/^j + J/^ PROCESS 

In this section, we calculate the branching ratio of the process r^^ — )■ J /ip + J /ip in the 
framework of NRQCD at one-loop level and in non-relativistic limit. Hence, the relations 
Pb = Pb = ^ : Pci = Pci = ^2^ and Pc2 = Pc-2, = ''2'^ are adopted. The bi-spinor operators 
are projected to states with the same quantum numbers as 775 and r^^, respectively, like 

and 

where iVc = 3, and Ic stands for the unit color matrix. In above, M^^ = Irub and Mjj^ = 2mc 
are implicitly assumed. 
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FIG. 1: Typical Feynman diagrams of the exclusive process rjh{Prj^) — )• J/ip{Pjf^^) + J/tp{Pj/ip2) 
at the one-loop level. 



For this process, at the leading order of relative velocity v in the framework of NRQCD, 
the tree level feynman-diagram has no contribution to the branching ratio, since the trace of 
the b-quark line form a Lorentz antisymmetric tensor, while the trace of c-quark line form a 
Lorentz sjTiimetric tenser. This situation remains also in the NLO counterterm, self-energy 
and vertex correction diagrams. Therefore, at one-look level only a few types of Feynman 
diagrams should be taken into account in the calculation, which are schemetically showed 
in Figure [H 

Because of parity and Lorentz invariance, the decay amplitude possesses the following 
unique tensor structure: 



A^(Ai, A2) - Ae^i,p„e*fi^^{\i)eJi^^{\2)Pji^^Pj/^^ ■ 



(3) 



In our calculation, the computer algebra system MATHEMATICA is ernployed with the 
help of the packages, FYENCALC I22J, FYENART [23i], and LoopTools [24. FYENART 
is used to draw Feynman diagrams and generate amplitudes, FYENCALC is used to trace 
the 7 matrices and to reduce various Passarino-Veltman tensor integrals 25|] to scalar ones, 
LoopTools is used to evaluate these integrals. For the aim of comparison, we express all the 
Passarino-Veltman scalar integrals encountered in this calculation in the Appendix. 

After taking the above mentioned procedures, it is straightforward to calculate this 
process and obtain the analytical amplitude in reduced form, i.e., 
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where 
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--Do(2)m^ + 2Do(l)m>,2 + Do(2)m>2 + —Co{l)ml + -C^{2)ml 
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Here, the form factor F (m^, m^) is a complex function; Dq, Cq, and Bq represent four-point, 
three-point and two-point Passarino-Veltman scalar integrals, respectively. The real part of 
F {ml, ml) is too complicated to be presented here, and therefore only the asymptotic form 
in small m,c limit is given: 

in 1 c: r 

Re(F {mlml))asy = ^log'(a) - -log(2)log(a) + -log(a) + ^log'(2) 

1, ,„, 297r2 3^3 3 

+ -log 2 + ^—n + - (6) 

2 ^^ ^ 96 8 4 ^ ^ 
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with a = ^. The full imaginary part of F {m,^, m^) is 

h 

^ ,„( 1 2^^ (36a-19)7r ,1 + 5^ (36a-5)7r, ,3 -(5, 



where 5 = a/1 — 4a, and its asymptotic form in the small nic limit reads 

XQtj- 77J- 

Im(F (m2,m^)),,, = — log(a) + Y^log(2) + vr . (8) 

With the above preparation, we can readily obtain the branching fraction of the exclusive 

"Qh -^ J I'^J 1'^ decay process. 



Br^.^mJliA-K,i ^^^^^^>*^^ Fi^,^. (9) 

Here, the dominant r^^ gluonic decay width is taken to be its total width approximately at 



one-loop order 



I.e., 



iTxal 2 



n^Atotai ~ r^z.o[r76 ^ gg\ = ir,,^^4(o) (lo) 

with 

r^ T , /^ / p. , TT^ 199 137r2 16 .as{2mb) 

Kgg = l + {Cf{-5 + ^) + Ca{— - -^) - -^nfTp) . (11) 



In numerical calculation, the following inputs are adopted: 

^Pj/^{0) = 0.263 GeV^/^ m^ = 1.5 GeV, mf, = 4.7 GeV, a, = 0.18 ~ 0.26 , (12) 

where the radial wave function at the origin ?/'j/^(0) is obtained by fitting the NLO QCD 
calculation result to the J/ip di-lepton decay width [27]. With the above preparation, we 
can readily obtain the numerical result of the concerned process 

Br[rit, -^ J/^J/tp] = 5.93 x 10"*^ ~ 2.58 x 10"^ . (13) 

Here, the uncertainties are originated from energy scale variation from rrib to rric. It is worth 
emphasizing that in our result the double logarithms exist and agree with what obtained 
in Ref. |8|, whereas our constant term does not agree with theirs, though its numerical 
influence is no big. 

III. HIGHER TWIST CONTRIBUTIONS 

As mentioned in preceding sections, at Born level r(?7j, -^ J/ipJ/ip) is exactly zero 
in non-relativistic limit, while the NLO radiative corrections are very small. People find 
that although in light cone formalism the leading twist term in the light cone distribution 
amplitudes(LCDAs) for J /ip vanishes in rjb — )■ J jipj /ip process, contributions from higher 



twist terms seem to be important [9]. In Ref. [9] the LCDAs up to twist-4 are taken into 
account for the consistency reason. It is true and in the following we reevaluate this process 
also in the light cone framework. However, to execute the twist expansion accurately, for final 
vector mesons with transverse polarizations, we expand the LCDA projector in momentum 
space given by 28| to twist-4, which yields more terms than what employed in Ref. J9|. i.e.. 

Ml = {Mf^ + Mf^ + M^) , (14) 
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Mf = ^f^E^^^^M^) , (15) 
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with E = [p^ + |p|)/2 and the transverse polarization vector 
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The higher twist LCDAs are related to the twist-2 ones by the Wandzura-Wilczek 
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iere, the contributions from three-particle DAs have been neglected as performed in Ref. 
9|. 

After a lengthy calculation, we get the final expression for the decay amplitude of the 
process rjQ -)■ V1V2, 

1 



M±± = Toef,^p^e/^el''^n''_n1 / dui / du2 
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Here, iprjoi^) is the wave function at the origin for pseudoscalar tjq. Our analytical result is 



I 



different from what given in Ref.[9|, partly due to the different projectors used. 
With the asymptotic form for twist-2 distribution amplitudes, 

<P±{u) = 0||(m) = (t)Asiu) = 6m(1 - u) , (27) 

the analytical decay amplitude turns to be pretty simple, it reads 

(vr^ - 4)mv,mv2{fvJv2 + fvjvi) 
+ ^A<JvJi+<fvJv.) ■ (28) 



-M±± = ne,^,^el'ie;in^_nl^^^^^ x 



To numerically evaluate the branching ratio of r]b — t- J/tpJ/ip process, we use the following 
commonly accepted input parameters: the charm quark mass in the MS scheme, m^^^ = 1.2 
GeV; the J/tb decay constant fj/^ = 416 MeV; and the /J/ . is obtained in the framework 
of NRQCD [29], /J/^ = 379 MeV. Then the numerical result reads 

Br[rib -^ J/ipJ/i^] = (1.1 ~ 2.3) x 10^^ . (29) 

Here, the uncertainties are also induced by the scale variation from mi, to rric as in above 
NLO QCD calculation. Note that the above magnitude is bigger than what the NLO result. 
This is however understandable, since roughly speaking the higher twist contributions are 
suppressed by factor of ( ^ j while the NLO contributions are suppressed by a^ ^ j . 

Following we apply the above higher twist analysis to the rjc to light vector mesons decay 
process for the first time to twist-4. As mentioned in the introduction the disagreement 
of experimental measurement with theoretical expectation is a long lasting issue. Before 
attributing some non-perturbative scenarios, to evaluate these processes in light cone for- 
malism to twist-4 is still meaningful. We already know that the leading twist term in LCDAs 
at leading order of as does not contribute to these processes, and the radiative corrections 
are dramatically suppressed by the light quark masses. Therefore, it is obvious that contri- 
butions from higher twist DAs dominate over others in the framework of perturbative QCD. 
In addition to the asymptotic form, the LCDA form in terms of Gegenbauer polynomials is 
also employed. That is 

0||,x(w, /i') = Qu{l -u)ll + Y^ a|l'^(/i2)C^/2(2M - 1) j . (30) 



The input parameters needed in the numerical calculation are listed in TABLE. HI fyS come 



from the QCD sum rules 30l-l33| . and the reasonable values of Gegenbauer moments ai and 02 
are from Ref . [3J] . All scale-dependent quantities refer to fj, = 1 GeV. The numerical results 

TABLE L Summary of theoretical input parameters. 

p K* LO (j) 



my [MeV] 


770 


892 


782 


1020 


/v[MeV] 


205 ±9 


217 ±5 


195 ±3 


231 ±4 


/^[MeV] 


160 ± 10 


170 ± 10 


145 ± 10 


200 ± 10 


ll>-L 
«1 





0.10 ±0.07 








02 


o.o9t°:J? 


n7+o-09 


o.o9l°:S? 


0.06+°:°? 



are given in Table. [TTl where Br[AS] and Br[GP] represent the results for forms of asymptotic 
and the Gegenbauer polynomials in LCDAs, respectively. Since the Gegenbauer moments 
a„ are small, the numerical results are not sensitive to the form of the leading distribution 
amplitudes, which influence the higher twist results via relations (iT9ll -( l2T]) . From results in 
Table [TTl we see that although the higher twist effect are tremendous for Bi[ric — )■ VV], it is 
still not enough to explain the experimental data. 

TABLE II: Experimental data and Numerical results for Br[r]c — )• VV], experimental data are from 
Particle Data Book [23]. 



Final state 


Br [ex] 


Br [AS] 


Br[GP] 


PP 


(2.0 ±0.7) X 10-2 


2.0 X 10-^ 


2.8 X 10-"^ 


K*K* 


(9.2 ±3.4) X 10-3 


7.2 X 10-^ 


9.0 X lO"'' 


UJUJ 


< 3.1 X 10-3 


9.1 X 10-5 


1.3 X lO-'' 


<j)4> 


(2.7 ±0.9) X 10-3 


6.6 X 10-^ 


8.1 X 10-4 



IV. CONCLUSIONS 

To further study the nature of recently observed state rjb, direct measurement of its 
decay products is necessary. The 77^ — )■ J jipj /tp process was considered and suggested to 
be a golden channel to this aim. In the literature, different theoretical estimation varies 



greatly. The branching ratio starts from 10~^ to 10~*, which induces some confusion for 
future experimental test. Within the pQCD and factorization scheme we have calculated 
this helicity conservation suppressed process, the rjf, — )■ J/ipJ/ip, at the next-to-leading order 
in pQCD. Our result confirms the existence of double logarithms, the log^(^) in Ref. 



-b 

^2 , 



and the coefficients of both double logarithm log^(-^) and single logarithm log(^) in our 

b b 

calculation are consistent with those in the same reference. However, we find that other 



B 



terms in our result deviate from those in Ref. |8| , though the numerical significance of the 
difference is not high. 

In the light cone formalism, the leading twist contribution to r]h — )■ J/ipJ/tp process 
vanishes. In this work we also evaluate the higher twist contributions to it. Expanding the 
LCDAs of final vector mesons to twist-4, we find that the higher twist terms contribute more 
to the decay width than what from the NLO corrections, which implies that the final state 
mass effects is more significant than the NLO corrections in this helicity suppressed process. 
According to our twist-4 calculation, the branching fraction of r^^ -^ J/ipJ/ip process can 



be as large as ^ 
feasible. In Ref. 



0~ , which enables the direct search of ?7{, in Tevatron Run II or LHC 
91], the same process was evaluated in the light cone formalism also to 
twist-4, but with different twist expansion procedure, which lead to different LCDAs from 
ours. We believe what we used are generated from the LCDA definition in twist expansion 
and should be more proper. 

Unlike the undiscovered rjf, — )■ J/ipJ/ip process, experimental results about rjc — )■ VV 
indicate that relatively large violations of the helicity conservation rules exist in these pro- 
cesses. The surprisingly large branching ratios of rjc ~^ VV still stand as a bewildering 
puzzle in charmonium physics. Since it is believed that the NLO corrections are greatly 
suppressed by the light quark mass, the higher twist effects might be large. In the light 
cone formalism, we have calculated this process with taking the next-to-next-leading twist 
effects in LCDAs of final vector mesons into account. Result shows that the higher twist 
DAs indeed violate the helicity conservation rules, though it still deviates a lot from the 
experimental measurement. This implies that the perturbative description of r^^ decay alone 
is not enough, and some non-perturbative mechanism m ay p lay important roles in 7]c decays, 



such as r]c — )■ VV, which deserves further investigation 35 1. 
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Appendix 

In this appendix, we list various Passarino-Veltman scalar integrals appearing in Eq.dS]), 
and only the leading power in a for real part are extracted; while for imaginary part, we 
present the full expression. Since coefficients of the Do(l) and Co(4)/m^ in Eq.(|5]) are same 
up to a sign, hence we only present their difference. Our evaluation for these integrals agrees 
with Loop Tools, and we have also checked with FIESTA2 [36| the asymptotic expression. 
Here, a = ^ and 5 = \/l — 4a. 



^' 



Do(l) = Do[m^, m^, 4m^, m^, m^, 2m^ + m^, 0, 0, m^, m; 
_ Co(4) log(2) 
ml mlml 



(31) 



Do(2) = DoK, ml, Ami ^l, "^^l + ^l^ '^^l + ^l, ^l^ ^l^ 0, 0] 
1 — log(a) .1 



Co(l) = Co[mj,, m^, m^, 0, 0, m^ 



31og2(a)+vr2 . log(i±f^ 



6m^ ml5 



(33) 



Co(2) =Co[m2,m2,2m2 + m2,0,0,m,2] 



-61og(2)log(a) + 31og2(2) + 7r2 . log(|^ 



2(o\ I ^2 Ino-rSz:^^ 



Qml ml6 



i^^^T^ (34) 



Co(3) =Co[4mf,m^,2mf + m^,0,0,mf] 



-61og(2)log(a) + 31og^(2) + 7r^ . log(fej: 

K ITT TTT- 

12ml 2m^(5 
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(35) 



Co(4) = Co[m^, m^, Am^,m^, 0, m, 



(36) 



Co(5) =Co[m^,m^,m^,m^,m^,m^] 



vr^ . log(2-4a) 

— ITT- 



12ml 



mi5 



(37) 



Co(6) 



Co[ml,ml,2ml + ml, ml, ml, 0] 
-61og(2)log(a)-31o g^(2) + 7r^ . log(^^^§^ 

ITT- 



6m^ 



m^(5 



Co(7) = Co[m5,m^,m^,m;,,m{,,mJ 



TT 



9m^ 



Bo(l) =Bo[m2,0,0]^- + 2-log(m2) + i7r 



Bo(2) = Bo[4m2, 0, 0] = - + 2 - log(4m2) + ivr 



Bo(3) = Bo[2mj^ + ml, 0, m^] ^ - + 2 - log(2m2) + Itt 



Bo (4) = Bo [m^,mj,,m^ 



1 - loEcfm. 



TT 



Bo(5) = BoK, m^, ^2| ^ _ + 2 - -)= - log(m^) 



Bo(6) =Bo[ml,0,ml] = - + 2-\ogiml 



2 + a 



(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
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